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Abstract 

Minkowski functionals constitute a family of order parameters which discriminate spatial 
patterns according to size, shape and connectivity. Here we point out, that these scalar de- 
scriptors can be complemented by vector-valued curvature measures also known as Quermafi 
vectors. Using examples of galaxy clusters, we demonstrate that the Quermafi vectors pro- 
vide additional morphological information on directional features and symmetries displayed 
by spatial data. 

Physica A 

PACS: 02. 40. Ft (conv. sets and geometrical inequalities), 02.40.Vh (global analysis and 
analysis of manifolds), 98.65. Cw (clusters of galaxies) 
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Figure 1: An isodensity contour of the galaxy distribution within the PSCz survey (from the 
PSCz homepage at http:/ /star-www.dur.ac.uk/cosmology/theory/pscz.htm] ). 



1 Introduction 

Spatial patterns originating from the polymorphic aggregation of matter in Nature occur on 
vastly different length-scales and with unlimited variety. The pattern shown in Figure |l| displays 
an isodensity contour of the smoothed galaxy distribution in our Universe, but it could also 
represent a biological structure obtained with X-ray tomography, for instance. 
The visual impression of such patterns may already convey valuable insights on the formation 
processes or the functionality of those structures. However, for unbiased inferences as well as for 
comparison with model simulations, for example, one needs objective and quantitative measures 
to characterize the geometry and topology of typical structural motifs. 

In the present paper we discuss a family of suitable measures, which include, besides the scalar 
Minkowski functionals, their vectorial relatives - the curvature centroids, or Quermafi vectors. 
Taken together, these measures constitute a complete class of morphometric order-parameters 
for spatial patterns formed - or approximated - by the set union and intersection of an arbitrary 
but finite number of convex "pixels" or "voxels". 

The usefulness of Minkowski functionals has already been demonstrated by the analysis of galaxy 
catalogues and simulations of large-scale structure [11, 1^, 16, [18[|, as well as of the temperature 



distribution in maps of the microwave background radiation |19[]. Furthermore, they were used 
to characterize the morphology of spinodal decomposition |15] and of patterns arising within 



reaction-diffusion systems |13]. Here, we focus attention on the less familiar Quermafi vectors 
and illustrate their versatility with the examples of galaxy clusters and with patterns generated 
by lattice automata. 
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2 Morphometric order parameters 



In this section we recall the definitions of Minkowski functionals and QuermaB vectors. The 
presentation will be brief and informal. Readers, who are interested in mathematical details, 
may consult refs. |, |, [l| . 



Consider a compact convex point set {body) C E'^, in Euclidean d-dimensional space with 
volume 



K 



VoiK)= I d'^x>0 (1) 
and with center of mass at the position: 

PO = Vo/Vo , Vo{K) = I xd'^x . (2) 



K 



Now, we let the body grow to form a parallel body K^, e > 0, by including all points x within a 
distance d{K,x.) < e to K. From Steiner's theorem ||l2| we know that the volume Vq^K,:) of iiTg 
is a finite polynomial in e. 



Vo{K,) = Y,(-]^MK)e^ > (3) 

i=o 

where = 7r*/^/r(l + i/2) denotes the volume of the i~dimensional unit ball. The coefficients 
Vi{K) define the Minkowski functionals of the body K. Since the volume Vq{K^) depends on 
the shape of K, the functionals Vi{K) contain morphological information on the original body. 
A corresponding expansion holds for the vectorial 'Vq{K) 

YoiK,)=Y,(%^^^iKy , (4) 

i=0 

which defines the Quermafi vectors \'i{K), < i < d. Under translations (t) and rotations (r) 
of the body K, the Minkowski functionals are motion-invariant scalars, whereas the Quermafi 
vectors are motion-equivariant: 

Yi{tK)= Yi{K) + tV^{K), (5a) 
V,(rK) = R{r)\,{K), (5b) 

where R{r) denotes an orthogonal rotation matrix and t a translation vector. Both the scalars Vi 
and the vectors Vj are additive functionals in the following sense: If the point set constituting 
the body K d'E/^ is the set-union of two bodies Ki and K2, then 

^i{KiyjK2) = 0i{Ki) + 0i{K2)-0i{Ki^K2) , (6) 

where we use the compact notation (p = {V^, Vj}. If Ki and K2 are arbitrary compact convex 
bodies, then the union set Ki U K2 is generally no longer convex. However, the right-hand- 
side of Equation (^) is still well-defined and may be taken as the definition of 0i{Ki U ^^2)- 
Consequently, additivity allows us to extend the application of the functionals to patterns 
formed by an arbitrary finite union of convex bodies via iteration of Equation (^). 
The functionals 0i, for 1 < i < d, depend on the shape of the body K; if its surface, dK, is 
smooth, with principal curvatures Cj(x), 1 < j < d — 1, then these functionals may be expressed 
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Figure 2: Examples of centroids for a series of simple patterns. The panels (l)-(4) show how the 
curvature centroids reflect the symmetry of a "cluster" . For patterns consisting of more than 
one body, the centroids weight the components according to their scalar Minkowski functionals. 
Note that the patterns in panel (5) and (6) can be discriminated using the centroids, but not 
by merely employing the scalar Minkowski functionals. 



in terms of surface integrals. For instance, in d = 3, with the surface area element dS'(x), and 
z = (l,x) we have: 

0^{K) = l[ zd5(x) , 0,iK) = ^[ i(ci(x) + C2(x))zd5(x) , (7) 

O JdK -^TT JdK ^ 

(psiK) = — ci(x)c2(x)zd5(x) . 
4vr JdK 

The vector V3(i^) is known as the Steiner point of the body K. - We list the meanings of the 
Minkowski functionals and the corresponding vectors in two dimensions in Table |l[ 
The most prominent feature of the family {(^j(-)|0 < i < d} is expressed by Hadwiger's char- 
acterization theorem: Consider the polyconvex ring P*^ of spatial patterns, defined to comprise 
the collection of arbitrary finite unions of convex bodies, /Cat = U^=i G P"^) N < oo, 
< diuiKa < d. According to the theorem, any additive, equi- or invariant and conditionally 
continuous^ functional $ on P'^ is uniquely determined by the linear combination: 

d 

HICN)=J2a,0,{ICN), (8) 

with the real coefficients Oj independent of the pattern /Cat. In this sense, the family {<fi} forms 
a complete basis of morphometric descriptors. 

For the practical analysis of cluster structures, for instance, it is convenient to employ the 
(curvature) centroids Pi = {pf,p^, ■■■) = ^i/Vi, provided Vi ^ 0, i = 0, .., d. 

^The conditional continuity requires that a functional is continuous with respect to the Blaschke-Haussdorff 
metric on the subset K'* C P'* of convex bodies. 
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Figure 3: Two initial configurations for the "Game of Life" (middle panel) and the results for 
the 188th generation (right and left panel). In the initial state, the "thunderbirdfuse" (black 
cells) is perturbed with the grey cell marked by an arrow. Whereas the "thunderbirdfuse" ends 
with a "blinker" (left panel), the "perturbed thunderbirdfuse" evolves into a pattern ejecting 
two gliders (right panel). 
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Figure 4: The Minkowski functionals for the two "Game of Life" series depicted in Figure ^. We 
compare the length of the circumference (left panel) and both components of the corresponding 
centroid (middle and right panel), "thunderbirdfuse": solid line, "perturbed thunderbirdfuse": 
dashed line. 



Let us consider the series of 2d-clusters shown in Figure ^ to gain some intuition on the kind of 
information supplied by the functionals (pi, and the centroids, in particular. 
The clusters (l)-(4) display the effect of symmetry reduction. In the case of the disk (1) and the 
ellipse (2), the centroids coincide at the symmetry center. However, Ci may be distinguished 
from C2 by the value of the isoperimetric ratio, I{Ci) = 1 (which holds only for the disk) and 
/(C2) > 1, where I{K) = 4V^{K)/[7rVo{K)V2{K)] > 1 for convex bodies. The centroids of C3 
and (^4 no longer coincide; in C3 they fan out along the single mirror reflection line, whereas 
in C4 they form a triangle. The clusters C5 and Cq have the same three components, but in 
different configurations. The scalar Minkowski functionals of both C5 and Cq differ from the 
ones of C1-C4 at least by the value of the Euler characteristic V2, which counts here the number 
of cluster components. But the scalar Minkowski functionals are unsuitable to distinguish C5 
from Cq - additivity implies ^^^(Cs) = V^Cq) ~, whereas the centroids discriminate clearly. 
In Figure ^ we illustrate the order parameters using a dynamical cellular automaton generated 
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Figure 5: The interplay of all curvature centroids for the second "Game of Life" series depicted 
in Figure ^ allows us to answer questions on the symmetry and movement of the patterns arising. 



Minkowski functional 


meaning 


corresponding vector 


Vo 


surface content 




Vi 


one fourth of the length of perimeter 


4 d 5'x 


V2 


Euler characteristic 





Table 1: The meaning of the Minkowski functionals in two dimensions, c(x) is the local curvature. 



by the "Game of Life" ||5|, |^. We consider two series of patterns, the "thunderbirdfuse" Q and a 
series starting with a slightly perturbed initial configuration of the "thunderbirdfuse" , with one 
point added (left column of Figure ^) . The "thunderbirdfuse" lights consecutively the "blinkers" 
visible as bars. The temporal sequence as well as the spatial course of the ignitions is reflected 
by the Minkowski functionals (Figure and the curvature centroids. The evolution terminates 
with "traffic lights" oscillating with period two and constant values of the scalar and vectorial 
functionals. 

The evolution of the slightly perturbed "thunderbirdfuse" proceeds rather differently. The early 
patterns are still comparable; at the time, when the pattern reaches the modified pixel, the 
perturbed "thunderbirdfuse" grows and attains a constant surface content from generation 188 
on. However, as seen from pi, the pattern keeps changing by ejecting two gliders visible in 
the lower right panel. The behaviour of the curvature centroids is shown in Figure ^: The 
asymmetry of the pattern and its temporal variation (the "gliders" move one cell in diagonal 
direction within four generations) is clearly recognizable. 



3 Applications to galaxy clusters 

Evidently, observational data of galaxies within clusters or X-ray photon maps do not provide 
from the outset the kind of spatial patterns required for the application of our morphometric 
tools. To start with, these patterns must first be constructed from the observational data; but 
there is no canonical way to proceed. Here, we employ two procedures, the excursion set method 



available, e.g. from the library of "Game of Life" patterns "xlife" (littp;//www. mindspring.com/~alanh/ 
life/). 
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Figure 6: Example for the excursion set approach. The data points of the cluster rCDM 
(middle panel) are smoothed using a smoothing length of 0.2/i~^Mpc. In the right panel we 
show the contours and the curvature centroids p, of the excursion set at a density threshold of 
u = 0.275N/ R'^. Here, denotes the number of particles within the cluster, Rm quantifies the 
cluster scale via the maximum distance of cluster particles from the cluster center. ~ h accounts 
for the uncertainty in the determination of the Hubble constant H: H = WOh^^^. - The left 
panel shows the counterpart of rCDM, OCDM, which is embedded into an open cosmological 
model. 



and the Boolean grain method in order to further illustrate our approach on the basis of simulated 
and real cluster data. 



3.1 The excursion set method 

We smooth the projected galaxy positions or the pixels of the X-ray photon maps with a Gaussian 
kernel. The smoothing length determines the scale of interest. Then we construct the excursion 
sets and investigate their topology and geometry using the Quermafi vectors and the Minkowski 



functionals [17| 



We illustrate this procedure by comparing two simulated clusters (part of the GIF-simulations, 
cf. 0]); they start with comparable initial conditions and evolve within different Friedmann- 
Lemaitre models as cosmological background, thus exemplifying the imprint of the background. 
We demonstrate our method in Figure ^ using the cluster rCDM from an Einstein-de Sitter 
model and a smoothing length of 0.2/i~^Mpc. The results of the Minkowski analysis are 
depicted in Figure |^. We plot the scalar Minkowski functionals vs. the density threshold 
defining the excursion sets: an averaged and smoothed density profile is encoded in the first 
Minkowski functional Vq; a comparison of the square of Vi and Vq quantifies how crooked the 
isodensity contours are. The Euler characteristic V2 counts the components (minus the number 
of eventual holes) of this substructure-rich cluster. 

The components pf of the curvature centroids are shown in the fourth panel of Figure ^. They 
wander in space if the density threshold is varied. This indicates that morphological features 
are shifting. Another way to make the centroids more illustrative is to consider the triangle 
spanned by the centroids and to compute its volume and the length of its perimeter (fifth panel 
of Figure |^). These quantities tell us how symmetric the isodensity contours are. 
To compare this cluster to OCDM, its counterpart in a low-density background model, we 

rCDM model is a variety of a Cold Dark Matter structure formation scenario embedded into an Einstein-de 
Sitter background. 
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Figure 7: First row: the scalar Minkowski functionals of the cluster depicted in Figure]^ vs. the 
density threshold u {u is given in units of the mean density within the cluster, is the number 
of cluster particles, Rm the maximum distance of cluster particles from the center of mass of all 
cluster particles). Second row: one component of the curvature centroids (left panel) and the 
area of the triangle spanned by the curvature centroids (right panel). 
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1^ o T* o Tvi ax ov 
pdi exilic Lei 


definition 




\J\jUL\L 


0.2 


C 


V((x-i)^)u 


3.20 


2.54 




(4yo(A(p))),/i?2 


8.78 X 10-^ 


6.31 X 10-"^ 


shifto /i?d 


rUmax dpi 




3.85 X 10"^ 


4.54 X 10"^ 



Table 2: The definitions of tlie condensed descriptors and tlieir values for the clusters rCDM 
and OCDM depicted in Figure 6. 4Vo(A(p)) is the circumference of the triangle spanned by the 
curvature centroids. Rd is the scale of the cluster, determined via the area of the convex hull of 
the cluster A^: ttR'^i = A^. a denotes the smoothing length. 



condense the morphometric information present in the Minkowski functionals and the centroids 
into a few dimensionless descriptors, by taking into account, that both clusters have different 
scales and numbers of particles. We define for a function / an average over density thresholds 
via: 



1 



if)u = —— I dn/(n) (9) 

"'max "'mi 



and consider the dumpiness, the symmetry parameter and the shift of morphological properties 
as defined in the third column of Table § 0. These parameters display the imprint of subclumps, 
the symmetry and shift of the curvature centroid p2 and are the larger the more substructure 
the cluster exhibits. As visible from Table ||, the cluster rCDM owns more substructure with 
respect to dumpiness and symmetry. A statistical comparison of these descriptors for larger 
numbers of clusters simulated in different models yields a typical cluster morphology and may 
help to constrain the present values of the cosmological parameters P]. 



3.2 The Boolean grain method 

The Boolean grain method decorates each point Xj of a point set with a ball i?Xi {''') of radius 
r. The union set of these balls is diagnosed by computing the Minkowski functionals and cen- 
troids as functions of the radius. To probe the morphology of an individual cluster locally, we 
place a window T) on the cluster center (the center of the point set) and study the centroids 
of By^.[r) n v. Inflation of the sampling window T) allows us to explore different regions of 
the cluster. We illustrate this method in Figure ^. Note, that also galaxies outside the window 
contribute ^ 

To provide a concrete example, we investigate the cluster CI 0016+161 at z ~ 0.54 observed by 
Belloni and Roser [Q] with the angular positions and spectral properties of galaxies in a field of 
3.5' X 5' and brighter than R = 23.5mag. Galaxies are considered as cluster members, if their 
redshifts lie within a predetermined rang^. We are left with 53 galaxies in this field, as shown 
in Figure |8|. 

We place spherical windows T) = 5/j^(c) of different sizes Rd on the center of mass c of the 
point set. The radius of the Boolean grains is held fixed at 15". In Figure |^ we show the 
behaviour of the scalar Minkowski functionals. The area of the union set of balls inside the 

''Here, Umin and Umax denote the maximum density found within cluster and outside the cluster, respectively, 
where "outside the cluster" means farther away from the cluster center than any cluster point 
^Computational details are described in 

®More precisely, following |^ elliptical and E+A galaxies are considered as cluster members, if for their redshift 
z: 0.525 < z < 0.575. For spiral and irregular galaxies the redshift range is 0.515 < z < 0.585. We consider only 
galaxies whose morphological type was determined (Table 5 in M]). 
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Figure 8: The cluster population of CI 0016+161 (left panel) and the associated Boolean grain 
model with circular windows (right panel). The radius of the Boolean grains is 15 arcseconds. 
The units are arcseconds. 




Figure 9: The scalar Minkowski functionals for the cluster CI 0016+161 vs. the radius Rd 
(in arcseconds) of the sampling window. The radius of the grains is fixed (15"), the Minkowski 
functionals are given in units of arcsec^"*. 
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Figure 10: The behaviour of the curvature centroids (pf-component) for the Boolean grain 
model in Figure ^ vs. the radius of the window (in units of arcseconds). The straight line 
indicates the cluster center. 



window V increases almost linearly with Rj:,. A homogenous galaxy distribution would yield a 
parabolic dependence, but here the galaxy density drops at larger scales. The small peaks in the 
Euler characteristic typically arise when two parts of a connected sub cluster enter the window 
simultaneously. 

The curvature centroids are more sensitive to geometrical anisotropics. The components are 
depicted in Figure IC. Two features are clearly visible: the centroids vary when the sampling 
window grows and, at a fixed scale of the window, significant differences between the curvature 
centroids arise. The variation of the centroids is clearly correlated with subclusters entering 
the window, e.g. the wandering of pi and po towards higher x-values reflects the large cluster 
component in positive x-direction of the center. 

To analyze the morphology of this cluster more in detail, we simulate a reference model for 
comparison. The simplest way to generate a reference cluster is an inhomogeneous, but spher- 
ically symmetric Poisson process We determine the projected density profile of the cluster 
non-parametrically around the center of mass using a binning, and simulate 100 Poisson clusters 
using this profile (the method for simulating inhomogeneous point processes is described e.g. 
in 1231). The mean Minkowski functionals and the centroids for the model clusters with their 



r.m.s. fiuctuations as well as the results for the real cluster are depicted in Figure 11 



The scalar Minkowski functionals deviate only slightly from the Poisson model. On the other 
hand, the curvature centroids reveal significant deviations (second row of Figure in par- 
ticular, the y-component of po (right panel) is shifted away from the true center of mass for 
relatively small values of the window scale, refiecting the fact, that there are more galaxies above 
than beneath the cluster center. 

- To strengthen our conclusions we also simulated the Poisson clusters with a different 
binning, but the results remained stable. The results demonstrate that the cluster C10016+161 
exhibits evidential subclustering in comparison with a Poissonian model. 
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Figure 11: The Minkowski functionals (first row) and components of different curvature cen- 

troids (second row) for the cluster CI 0016+161 (sohd hne) and for a set of Poissonian random 
clusters following the same density profile as CI 0016+161 (dashed line) vs. the radius of the 
Boolean grains in units of arcseconds. The radius of the window is fixed at 69". The shaded ar- 
eas indicate the one-cr-range for the random clusters. Whereas the scalar Minkowski functionals 
are consistent with the Poisson model, the curvature centroids, especially their j/-component, 
reveal a significant deviation from this model. 
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4 Conclusions 



The examples discussed in this note support our claim that the method of cluster analysis 
based on the combination of scalar Minkowski functionals with vectorial centroids furnishes a 
versatile set of order parameters to sort out essential aspects of cluster morphology such as 
symmetry, dumpiness, global shape and topology. Our cluster morphometry rests on a solid 
mathematical basis derived from a few reasonable requirements. In this sense, it allows for a 
unique morphological description. The construction of patterns from empirical data introduces 
additional parameters which may be employed advantageously for scale-specific diagnosis. No 
tacit statistical assumptions are involved. 

Finally, we note that these families of morphological measures may be extended further to 
include tensor-valued Minkowski functionals, which generalize the concept of inertia tensors. 
The code to compute the Minkowski functionals and the Quermafi vectors is available on request 
from the authors. 
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